Transmon-phonon coupling of plasma oscillations and lattice vibrations 
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In the transmon qubit we expect from conservation of momentum and energy a coupling between 
the plasma oscillations and the vibrations of the underlying lattice. Specifically, the electron veloc- 
ities and their kinetic energy density are boosted by the underlying lattice vibrations. We consider 
this effect in a representative transmon comprising two semi-circular superconducting charge islands 
joined by a Josephson junction. In particular, we solve the Fourier transform of a two-dimensional 
radial current density having inversion symmetry. The resulting spectral density is ohmic but also 
scales quadratically with the critical current I c and logarithmically with the size of the transmon: 
J(cj) ~ Ic ujlog(kR). We make positive-definite Born-Markov approximations in a generalized 
Fermi's Golden Rule and estimate the phonon-induced dephasing rate is negligible compared to 
current experiments. 



I. INTRODUCTION 

The transmon, or transmission-line shunted plasma os- 
cillation qubit, is a superconducting Josephson junction 
qubit situated in a microwave cavity^ - — The junction 
connects a pair of superconducting islands whose charg- 
ing energy Eq is significantly smaller than the Josephson 
energy Ej. The large ratio Ej/Ec ^S> 1 reduces sensi- 
tivity to charge noise relative to the anharmonicity and 
thus improves dephasing times relative to the speed of 
qubit operation. The use of a three-dimensional cavity^ 
and other refinements^ have improved dephasing times 
to almost 100 /j,s& 

The transmon is similar to the Cooper-pair box (CPB) 
qubit 7 £ but with a significantly larger ratio Ej/Ec- To 
achieve low charging energies the islands can be expan- 
sive, with dimensions measured in hundreds of microns 
and combined area approaching 0.5 mm 2 . It is typi- 
cally fabricated on a substrate such as sapphire, with 
the junction a thin oxidized layer such as AI2O3 between 
islands of superconducting Niobium or Aluminum. It is 
then placed in an electromagnetic cavity where it is mea- 
sured, controlled, and coupled to other transmons, by 
microwaves^ 

In the Ej/Ec 3> 1 regime the transmon is similar 
to an anharmonic oscillator (cosine potential) but with 
the number of excess electron pairs n and the supercon- 
ducting phase difference (f> as the conjugate variables. A 
superposition of the lowest eigenstates |0) and |1) will 
result in oscillations of the current I c sin(</>) in which the 
total momentum of the charge carriers is also oscillating. 
Conservation of momentum suggests a recoil of the un- 
derlying lattice. The relatively large area of the transmon 
also encompasses many lattice sites whose vibrations may 
interfere with qubit operation. 

In the case of superconducting quantum interference 
device (SQUID) qubits, the coupling to phonons was 
previously derived from the boost that lattice vibrations 
give to the electron velocities and thus their kinetic en- 
ergy density^ Here we apply the same technique to a 
representative transmon qubit using a radial model for 




FIG. 1. Model transmon comprising two semi-circular su- 
perconducting islands joined by a Josephson junction (bold 
x). Radial current density (thick arrows) flows off one is- 
land, through the junction, and onto the other. A phonon's 
transverse in-plane motion (thin arrows) collectively boosts 
the electrons' kinetic energy, especially within a diametrical 
band (dashed borders) aligned with the transverse motion. 



the current density. We solve its Fourier transform and 
find the spectral density is ohmic, but also quadratic in 
the critical current, and scales logarithmically with the 
size, in phonon wavelengths, of the transmon: J{uj) ~ 
1% ojlog(kR). We also use positive-definite Born-Markov 
approximations in a generalized Fermi's Golden Rule to 
estimate the phonon-induced dephasing time is orders of 
magnitude longer than state of the art experiments. 
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II. EIGENSTATES AND ENERGIES 

The total Hamiltonian comprises that of the transition, 
Hs, the phonon bath, Hg, and their coupling, V: 



H = H S + Hr 



V. 



(1) 



The Transmon CPB Hamiltonian 



Neglecting the microwave cavity, the transmon's effec- 
tive Hamilton can be reduced to CPB formi 



H s = AE c {n - n g ) 2 - Ej cos < 



(2) 



with <fr as the phase difference across the junction, n — 
—i-^p as the number of excess Cooper-pairs on one island, 
and n g as the effective offset charge in pair-charge units 
2e. 

By introducing the function g(x) = e~ 2mgX tp(2x) the 
eigenvalue equation becomes the Mathieu equation 



g"(x) 
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Ec 2Eq 



g(x) = 



(3) 



whose solutions me^^q, x) can then be used to reconstruct 
1 



= -=e ta 'W(j^/2). 
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The boundary condition ip(<fi) = ip(4'+2Tr) restricts the 
Mathieu functions me v {q 1 x) to characteristic exponents 
v = —2(n g — n c ) with n c an integer number of extra cycles 
of ip(4>) fitting into a full 2ir period of </> . The Mathieu 
characteristic eigenvalues a(u, q) = E/Ec thus separate 
into bands with n c = n c (m,n g ) sorting the eigenvalues 
according to band index m. 

As shown in Figure O the charge dispersions, or band 
widths, decrease exponentially with Ej/Eq while the 
transmon's relative anharmonicity, (E21 — Eiq)/E\q, de- 
creases only algebraically X Reducing the anharmonicity 
prolongs the time needed for qubit operations but the 
exponentially reduced dispersion more than compensates 
with extended coherence times. This is the advantage of 
the large Ej/Eq ratio in the transmon qubit. 

Within the {|0) = |V>o) > |1) = IV'i)} logical subspace 
we have H s ~ a z = |0) (0| — |1) (1|. Since the current 
through the junction 7 c sin0 couples to the phonons we 
will soon also be using the operator = sin </>. We calcu- 
late its non-zero (off-diagonal) transition matrix elements 
by numerical integration: 



<0|- 



ID = 



d(j) %{(j>) sm{4>)il)x{(t>)- 



(5) 



For Ej/Ec = 49, as an example, (0|cr^|l) rj 0.414« 
which is not far from the l/y2a ~ 0.450 that one would 
get from a small angle approximation sin(</>) f=a <f> in a 
not-anharmonic oscillator with a — \J (Ej / Ec)/8- 
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FIG. 2. Transmon energy levels and probability densities in 
a cosine potential. Left plot: lowest four bands with charge 
dispersion and anharmonicity for Ej/Ec = 9,25,49. Right 
plot: probability densities for Ej/Ec = 49. 



B. The Harmonic Phonon Bath 

The harmonic crystal Hamiltonian is composed of 
phonon modes labelled by wavevector k and polarization 
index s. Each phonon contributes an energy ?jwj s : 



at or . 

fes ks 



(6) 



k,s 



Here at and ar are the phonon creation and annihila- 

ks ks r 

tion operators. We assume linear dispersions, uJi^(k) = 
ct \k\ and u>z{k) = cl \k\, for transverse and longitudinal 
polarizations e s (k). The lattice site speeds are 



t h ^ks 

2pV 



t £ 

ks 



-ik-r 



e s (k) (7) 



with p and V the mass density and volume of the lattice. 

With the direction of phonon propagation k specified 
by angles 9 and tf> measured, respectively, from z per- 
pendicular to the transmon and x transecting it, and 
choosing ii as the transverse in-plane polarization, the 
polarizations have {x,y, z} components as follows: 

ii = (k x z)/\k x z\ = { sm.<j) — cos^ } 
62 = 63 x ei — { cos <j> cos 9 sin cf> cos 9 —sin 9 } 
£3 = k = { cos (p sin 9 sin sin 9 cos 9 } 



C. The Transmon-Phonon Coupling 



The coupling between the transmon and the crystal 
lattice arises from the fact that the current is formed 
from the electronic band states in the reference frame 
co-moving with the lattice sites during the vibrations of 
the crystal^ 

In the lab frame the electron velocity v e must therefore 
include the speed u of the lattice sites: v e — j / (en e ) + u. 
Here j is the current density and e and n e are the electron 
charge and number density of the electrons. Their kinetic 
energy density n. 



e o ' 



* e | thus acquires a cross term and 
the transmon-phonon coupling V accounts for this extra 
energy density integrated over the transmon: 



V 



e 



r j • u. 



(8) 



Qualitatively, the charge carriers converge radially in- 
ward from one semi-circular island, flow through the 
junction, and diverge radially outward onto the other 
island. The junction current, J c cr^, flows at a radius r 
through an area nrT with T as the thickness of the is- 
lands. Our current density model is thus radial above 
the cc-axis and opposite below and decreases like 1/nrT 
with distance r from the junction: 



Iccr<j, y 
irrT \y\ 



(9) 



Integrating j ■ u over the transmon amounts to solving 
the components jV of the Fourier transform of the cur- 
rent density. In the thin-transmon approximation (thick- 
ness T <C X/2tt for relevant phonon wavelengths A ~ 1 
//m) we need only integrate in the xy-plane of the trans- 
mon out to its radius R. In polar coordinates (r, ip), with 
ip measured from the a;-axis transecting the transmon, we 
define a dimensionless version of this Fourier transform, 



kit 



d(kr) I #-^f • e s (k)e^ krk - f (10) 



to write the coupling coefficients as 



9ks 



■Iks 
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(11) 



in terms of which the transmon-phonon coupling is 

v = ^ ® E {dksHs + ffL°L) ■ (12) 



III. THE FOURIER TRANSFORM 
A. Transverse in-plane polarization 

Transverse in-plane phonons, as shown for example in 
Figure [TJ are the relevant polarization for boosting the 
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FIG. 3. Directional emission pattern Ijgjl 2 of transverse in- 
plane phonons at kR — 5000. 



transmon's electron velocities. Whatever their direction 
of propagation, there is always a diametrical band within 
which the current flow is qualitatively aligned with the 
phonon's transverse in-plane collective motion. 

Quantitatively, we take advantage of the inversion 
symmetry to integrate cos(fcrfc • f)f ■ e% twice over one 
island. With k-f = cos(0 — tp) siri9 and r-e% = sxD.((j) — ip) 
we find 



Jki 



Si[fci? cos cj> sin 9] 



suit 



-4*5 (13) 



sin t 



where Si[ ] is the sine integral function. 

The corresponding directional emission pattern for 
transverse in-plane phonons, IjjJ 2 = 4 2 Si 2 [/cj;i?]/ sin 2 9, 
is shown in Figure [3] Steeply incident (0 <SC 1) phonons 
dominate the pattern because they have relatively long 
in-plane wavelengths A/sin# that give a wider band 
within which the coupling accumulates. On the other 
hand, phonons propagating in the ±y directions have 
their band of transverse motion straddling the islands' 
gap and the in- and out-flows cancel. Accordingly, the 
phonon-induced decoherence is mostly effected by steeply 
incident phonons in the x, z plane. 

Solving the spectral density will involve integration 
over all directions in fc-space. We obtain by numerical 
integration, for a variety of 250 < kR < 7500, a good fit 



2ji 



d'n = / d9 sin 

Jo 



with a w 496 and b w 1.16. 



lifcfwalogp. kR] 
(14) 
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FIG. 4. Ap arises from [V(ti), [Vfa), (caricature dark 

shading) integrated within the upright trapezoid. Symmetriz- 
ing the integral around £2 = ti discards self-adjoint shifts to 
Ho- Positive-definite Born-Markov approximations are made 
with a Gaussian cutoff in the relative time r and a commen- 
surate Gaussian smoothing in the average time T.— 
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IV. POSITIVE-DEFINITE BORN-MARKOV 
APPROXIMATIONS 

In this section we give an overview of techniques devel- 
oped by Taj et alfii to generalize Fermi's Golden Rule to 
nanodevices with positive-definite Born-Markov approx- 
imations. 

The unperturbed Hamiltonian Hq = H$ + Hb defines 
a standard interaction picture Ut = e~ lHot l h in which 
the coupling V(t) — UJVUt determines the evolution of 
the system-bath density operator p(t) = uJp(t)Ut as 

Ap = -- S / dti / dt 2 [V(h),[V(h),p(t2)}}, 

n Jt-At/2 Jo 

(15) 

where we have neglected a first order term whose contri- 
bution to the system density operator p$ will vanish, for 
a thermal bath, in the averaging A/55 = Trs[Ap] over 
bath eigenstates. 

Simply dividing Ap by At gives a course-grained ap- 
proximation to p ready for Born and Markov approxi- 
mations. However, the dynamics of p can be kept physi- 
cal (positive-definite) by an improved nested-integration 
around t\ = t\ — t w and t' 2 = ti — t ~ using a 
Gaussian cutoff in the relative time r = t\ — t' 2 and a 
Gaussian smoothing in the "macroscopic" average time 
T = + t' 2 ) /2 as followssii^ 



^1 ^2 

2 [u}V(T + 7/2)U t , [U}V(T^7/2)U t ,p(t)}} (16) 



Here we have made Born's approximation by replacing 
p(t 2 ) with pit) and Markov's approximation by using a 
Gaussian cutoff time f longer than the correlation times 
for the weak coupling V(t). In extending the relative time 
integral back to —00, we have also disposed of the self- 
adjoint parts of the double commutator structure that 
could just as well have been included as small shifts to 
the unperturbed Hamiltonian Hq. 

Finally, the smoothing time is chosen as T = f/2 so 
that the Gaussian product re-factors as 



-T 2 /2T 2 -r 2 /2f 2 -t' 2 /f 2 -t' 2 If 2 

e ' e ' = e 1 ' e 2 ' 



(17) 



in the t\ , t' 2 plane and the integrals separate in which case 
the evolution of the density operator may be written as 



P : 



2tt 1 
¥2 



VlLU u [U}LU u p]} 



in the interaction picture, or as 

2tt 1 



P 



h 



[Ho,p] 



h 2 2 



[L, [L,p\] 



(18) 



(19) 



in the Schrodinger picture, with the Lindblad operator 



L = Vw/V27T 3 / dt'V(t') 



(20) 



whose matrix elements between final and initial states are 
like the square root of a delta function as w e 1/r — >■ 0: 



Si</i^>i 2 = Si</m*>i 2 

h n 



-(u s -Ui) 2 /AQ 2 



V \/2ttuj 



(21) 



We thus see that Fermi's original Golden Rule is embed- 
ded in the more general positive-definite double commu- 
tator structure of [L, [L,p]]. 
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V. TRANSMON DEPHASING RATE Ti 

A. Master Equation 

The master equation for the transmon system is ob- 
tained by tracing over the bath, 



PS 



-jr[H s ,p s ] 



~Tr B [L,[L,p s ®PB}}, (22) 



where we start with the transmon ps initially uncorre- 
cted with a thermal bath p B = (e- HB / k ' BT )/Z B and 
extend Born's approximation by assuming the (small) 
changes in the interaction picture p occur mainly in the 
transmon: p(t) = ps(t)® p~B and thus p(t) = ps{t) ® p B - 

Recalling V = Y,%„ a <P ® (5fc s a £ s + 9^^) we solve in 
the two-dimensional logical eigenbasis of Hs the Lind- 

I 



blad operators L^ s for each mode of the bath, i.e. 



ks 







( 1 k*|0)5S.o s> o 



VS (23) 



where yS ~ ^/<5(oj^ s — oj), in terms of the transmon fre- 
quency cj = u)\ — ujq, and we have kept only the energy- 
conserving terms in the transitions between |0) and 
Taking L = J2k s L k s and writing 



PS <S> PB 



Poo Poi 

PlO Pll 



PB, 



(24) 



with p mn = (to | ps \n), we multiply out the 2x2 matri- 
ces and trace the bath, making frequent use of the cyclic 
property of the trace and Trs[at ap s ,ps] = 5^j:,5 ss ' N Uf . 
in terms of the thermal average occupation numbers 
= l/(e^/ fcBT - 1). We find 



Tr B [L,[L,p]] = \{0\, 



V 1 

(27r)3c7' 



d 2 a 



huj r 2 l 



e 2 2pV 7r 2 k< 



fell 



r 



2(N ulPoo - (N u + l)p u ) (2i\^ + l)p i 

(2A^ + l) Pw 2{{N„ + l)pn - N^poo) 



(25) 



where we have converted the fc-space sum to an integral, 
integrated over the delta function, and kept only the (s = 
1) in-plane polarization. Recalling the logarithmic fit to 
the numerical integration of the dimensionless Fourier 
transform, the spectral density reduces to 



2 rn 2 /e 2 
(2it) 5 ctP 



' huj a\og[b kR] (26) 



which is ohmic (oj is raised to the first power) but also de- 
pends quadratically on the critical current I c and grows 
logarithmically with the size of the transmon kR (mea- 
sured in phonon wavelengths). 



a corresponding model transmon radius of R = 400 pm. 
Operating at uj/2it = 4 GHz and Ej/E c = 49, with 
I c = 20 nA, and at low temperatures ksT <C fioj, e.g. 10 
mK, we calculate, with p ss 8570 kg/m 3 and ct ~ 1600 
m/s, a dephasing time of T^g 1 w 10 s. 

Operating a similar qubit 10 times as fast with 10 times 
the current, for example by using different junction pa- 
rameters and thus energy Ej, gives a dephasing time of 

7 ms which is still two orders of magnitude longer than 
IBM's almost 100 ps decoherence time. 



VI. CONCLUSION 



B. Thermalization and Dephasing 



We conclude that transmon phonons are no problem. 



In the master equation the (diagonal) population- 
mixing, or thermalization, vanishes in the thermal bal- 
ance 



Pn = 
Poo iV, 



_ e -hu/k B T 



(27) 



For the (off-diagonal) coherence poi we use (2iV w + l) = 
coih[fiLjj /2ksT] to write the dephasing rate as 

r 10 = ~\ <0| a |1> | 2 J>)coth[^]. (28) 

As an example, consider the IBM transmon^ whose 
niobium island/pads encompass sa 0.5 mm 2 . We choose 



Appendix: Other Polarizations 

For the other two polarizations we make use of a mul- 
tipole expansion of the periodic step function 



sin(ra^). 

7TJ7T 



(A.l) 



m— 1,3,... 

With f ■ {e 2 , £3} = cos(0 — ip){cos8, sin#} we find 



3k 



fe2,3 



y ± s inH)2« m - 1 #^^. (A.2) 



m— 1,3,... 



{tan 61,1} 
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We don't expect longitudinal phonons to have an ap- 
preciable effect over many wavelengths. The in-plane 
projection of the not-in-plane transverse polarization is 
similarly longitudinal. In fact, numerical integrations 
consistently show these spectral densities to be small 
fractions of that of the transverse in-plane polarization: 

J fei 2 » J <fn |.y 2 » J d?n |^ 3 | 2 . (A.3) 



We thus use only the transverse in-plane polarization in 
estimating the dephasing rate. 
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